Transport of ultracold Bose gases beyond the Gross-Pitaevskii description 
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We explore atom-laser-like transport processes of ultracold Bose-condensed atomic vapors in meso- 
scopic waveguide structures beyond the Gross-Pitaevskii mean-field theory. Based on a microscopic 
description of the transport process in the presence of a coherent source which models the outcou- 
pling from a reservoir of perfectly Bose-Einstein condensed atoms, we derive a system of coupled 
quantum evolution equations that describe the dynamics of a dilute condensed Bose gas in the frame- 
work of the Hartree-Fock-Bogoliubov approximation. We apply this method to study the transport 
of dilute Bose gases through an atomic quantum dot and through waveguides with disorder. Our 
numerical simulations reveal that the onset of an explictly time-dependent flow corresponds to the 
appearance of strong depletion of the condensate on the microscopic level and leads to a loss of 
global phase coherence. 
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I. INTRODUCTION 

The rapid progress in the experimental techniques for 
trapping and manipulating ultracold Bose gases on mi- 
croscopic scales has opened the possibility for investi- 
gating mesoscopic transport properties with interacting 
bosonic matter waves. A key achievement in this context 
is the development of simple or more complex waveguide 
geometries for cold atoms with optical techniques (as, 
e.g., in Ref. [l[) or with atom chips [1, H, 0], which also 
allow one to impose scattering and disorder potentials on 
microscopic scales (e.g. [HEl) and which also permit to 
detect individual atoms with rather good accuracy [in- 
coherent transport processes in such waveguides, where 
the atoms of the Bose gas are freely propagating along 
the guide with some hnite momentum, can be studied, on 
the one hand, through bosonic wave packets that are cre- 
ated upon sudden release of a Bose-Einstein condensate 
from a trap, which is then accelerated in thepresence 
of a finite potential gradient along the guide [a 0. On 
the other hand, the principle of an atom laser [T(| E| 
can be used for this purpose, as was recently demon- 
strated in Refs. [H, [lj, [3, 0j|. In these experiments, 
a coherent matter-wave beam was injected from a Bose- 
Einstein condensate in a trap into an optical waveguide. 
This was achieved by means of an radio-frequency (rf ) in- 
duced transition from a magnetically trapped (mp = — 1) 
to an untrapped (tof = 0) hyperfme state of 87 Rb in 
Ref. [H, EM, and by means of a careful ramping of a 
magnetic- field gradient in an optical trap in Ref. [lH ]. 
In this way it becomes possible to study bosonic scat- 
tering processes with cold atoms at well-defined incident 
energy, in close analogy to scattering of laser beams and 
to electronic transport in mesoscopic solid-state systems 

From the theoretical side, a number of investigations 



on the quasicontinuous transport of ultracold Bose gases 
and Bose-Einstein condensates have been undertaken 
during the past decade. This started with the attempt 
to define an atomic analog of Landauer's quantization 
of conductance and was continued by first inves- 
tigations of nonlinear resonant transport and interac- 
tion blockade in quantum-dot-like scattering potentials 
[3 ll^t - Propagation and transmission studies were 
then undertaken on the basis of the stationary Gross- 
Pitaevskii equation 0, HH, H3, [23[ as well as through a 
time-dependent integration approach in which the injec- 
tion of bosonic matter waves was accounted for by means 
of a coherent source term 0, [H, [26j] . 

With few exceptions (e.g., Ref. [19]), the above studies 
were mainly based on an effective one-dimensional Gross- 
Pitaevskii equation. This implicitly assumes the realiza- 
tion of a one-dimensional mean-field regime [27| where 
the transverse confinement of the waveguide is strong 
enough to inhibit the population of transversally excited 
modes, but not so strong as to come close to the Tonk- 
Girardeau regime of impenetrable atoms [HI, Hi| . Fur- 
thermore, the phase coherence length for a freely prop- 
agating atomic beam, which is generally finite for one- 
dimensional quasi-condensates [301 ] , is assumed to be con- 
siderably larger than the effective longitudinal extent of 
the waveguide (which may be limited due to a finite focal 
region in the case of an elongated dipole guide [12], or 
due to atom detection [7[). Using the one-dimensional 
Gross-Pitaevskii equation under these conditions, the 
presence of atom-atom interaction in the bosonic beam 
mainly manifests, as in nonlinear optics [3l| . in form fi- 
nite nonlinearity effects in the wave scattering process, 
such as nonlinearity-induced shifts of resonant transmis- 
sion peaks [ll[13,|2i]. 

The main focus of this work is to determine the quan- 
titative amount of depletion that is generated during the 
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transport process of the condensate, which is a rather rel- 
evant information from the experimental point of view. 
Indeed, too much depletion would eventually lead to a 
complete destruction of the phase coherence of the beam, 
which means that wave interference phenomena in the 
transport process, which would result from the coher- 
ent Gross-Pitaevskii equation, may, in this case, not 
be observable in practice. Inspired by previous stud- 
ies [HI, |H, [34[ , we use a quantum kinetic approach for 
this purpose, which is equivalent to the Hartree-Fock- 
Bogoliubov approximation [35| . In practice, the conden- 
sate wavefunction is, in this approach, propagated to- 
gether with two-component functions that describe the 
coherence as well as the one-body density matrix asso- 
ciated with non-condensed atoms. This ansatz generally 
involves a renormalization of the microscopic interaction 
strength, which, however, does not represent a concep- 
tual problem in the quasi-one-dimcnsional confinement 
geometries that we are considering here. Moreover, it ex- 
plicitly breaks the gauge symmetry of the Bose ga s and 
introduces a gap in the excitation spectrum [34], [35[ . We 
justify this due to the fact that the matter-wave beam 
is, in the guided atom-laser scenario, connected to an 
idealized reservoir that contains a macroscopically large 
Bose-Einstein condensate, for which spontaneous break- 
ing of gauge symmetry can be well assumed [36[ . 



In order to provide a solid theoretical foundation of 
our approach, we begin, in Section [TH with a step-by- 
step derivation of the quantum kinetic equations that are 
used to model the propagation process of the condensate 
in the presence of depletion. For the sake of definitcness, 
we focus here on an atom laser process that is based on 
the rf-induced outcoupling of atoms from a Bose-Einstein 
condensate in a magnetic trap [T2I Il3| . Our starting 
point is the microscopic many-body description of the 
two-component system that defines this atom laser, from 
which we derive, under the assumption of a weak effec- 
tive atom-atom interaction in the waveguide, nonlinear 
kinetic equations for the condensate wavefunction and 
for two-component functions that account, in lowest or- 
der, for quantum depletion. In Sections IIIII and IIV1 re- 
spectively, we apply this approach to transport processes 
of Bose-Einstein condensates through double barrier po- 
tentials and disorder potentials that were also studied in 
Refs. [U, [25| on the basis of the Gross-Pitaevskii equa- 
tion. Our main finding, which is summarized in the con- 
clusion in Section [V] is that the overall coherence of the 
atomic beam remains fairly well preserved over reason- 
ably long time scales in the case of a quasi-stationary 
flow of the condensate, while strong depletion arises in 
the case of permanently time-dependent scattering in the 
waveguide. 



II. KINETIC EQUATIONS FOR THE 
SCATTERING SYSTEM 

We consider a gas of ultracold bosonic atoms that are 
confined in a large trap and injected from there into 
a mesoscopic waveguide structure. In practice, such a 
"guided atom laser" can be realized by combining a mag- 
netic trapping potential with an optical waveguide which 
is, e.g., created by an elongated dipole trap. A radio- 
frequency (rf) field is then used to flip the spins of the 
atoms and induce coherent transitions from a low-field 
seeker state |r) (such as the \F, wif) = | X, — 1^ state of 
87 Rb) to a hyperfine state \g) that is, in lowest order, 
insensitive to the magnetic field (the \F,mp) = |0, 0) 
state of 87 Rb) but still experiences the optical potential. 
As demonstrated in Refs. [H, 03 where this setup was 
experimentally realized, the quadratic Zeeman effect can 
be exploited to compensate the quadratic variation of the 
laser intensity along the dipole trap and to thereby create 
a rather homogeneous waveguide potential. 

In the formalism of second quantization, the quantum 
dynamics of this two-component system is described by 
the coupled equations 

ih^ g (r,t) = (-^-A + V g (r)^ g (r,t) 

+tf rfl tft(r,i)tf r (r,i)tf fl (r,t) 
+tf(t)* r (r,t) (1) 

ih^r(v,t) = f-^A + K,(r)W(r,i) 
at \ 2m J 

+U r ¥ r {r,t)i> r {r,t)% r {r,t) 

+U rg ¥ g (r, t)4f g (r, t)* r (r,i) 

+K'(t)9 B (T,t) (2) 

for the bosonic field operators ^> g (r,t) and ty r (r,t) that 
annihilate atoms in the "waveguide state" \g) and the 
"reservoir state" |r), respectively. V g (r) denotes the 
waveguide potential while V r (r) corresponds to the mag- 
netic trap (perturbed by the presence of the optical 
guide) that confines the reservoir atoms in the state |r). 
The short-range interaction between the atoms in the in- 
volcd hyperfine states is represented by contact potentials 
of the form U{t\ — r%) cx 8{r\ — r^) with the prefactors U g , 
U r , and U rg for \g)-\g) 7 \r)-\r), and |g)-|r) interaction pro- 
cesses, respectively. Transitions between waveguide and 
reservoir states are induced by the coupling amplitude 
Kit), which is adiabatically raised from zero to a finite 
value and which would correspond to the field strength 
of the rf radiation. 

It is reasonable from the experimental point of view to 
consider an initial many-body state in the trap that con- 
sists of a nearly perfect Bose-Einstein condensate con- 
taining a rather large number of atoms. Provided the 
coupling amplitude K(t) is not too strong, we can safely 
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assume that the condensate is not appreciably affected by 
the outcoupling process on finite time scales. We there- 
fore make the ansatz 

$ r (r, t) = (# r (r, t)) = r (r) exp(-i/it/») (3) 

where ^(r) is the macroscopically populated condensate 
wavefunction that satisfies the Gross-Pitaevskii equation 



magnetic trap 



A + V r (r) + U r \y r (r)\ 2 ) * r (r) = //* r (r) (4) 



2m 

with the normalization 



d 3 r|* r (r)| 2 =jV P » 1, 



(5) 



and \i denotes the chemical potential of the condensate. 
This ansatz explicitly involves a spontaneous breaking of 
the condensate's gauge symmetry, which should be valid 
in the thermodynamic limit of a very large number Af r 
of atoms in the trap [361 ]. Inserting this expression for 
\& r (r,t) into Eq. |Q} yields then 



-A + V g (r) + U rg \y r (r)r)V g ir,t) 



2m 

U g fy{T,t)%(T,t)* g (T,t) 

if(f)* r (r) exp(~ifit/h) 



(6) 



as equation for the field operator of the waveguide state, 
which contains a coherent source term with the ampli- 
tude K(t)# r (r). 

The external potential for the atoms in the state \g) is 
written as 
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V g (v) = -muiri + V(x) 



(7) 



where x and rj_ = (y, z) denote the spatial coordinates 
along and perpendicular to the waveguide, respectively, 
and u± is the transverse confinement frequency. In addi- 
tion to the waveguide, we consider the presence of a lon- 
gitudinal scattering potential V(x) of finite spatial range, 
which could, e.g., be a sequence of barriers or speckle dis- 
order induced by other laser fields. Denoting the trans- 
verse eigenstates within the waveguide by x n (rj_), we can 
make the decomposition ^ g (r,t) = J2 n Xn{r±)?pn(x,t) 
for the field operator of the |<7)-atoms, where ip n (x,t) 
annihilates a particle at position x in the nth trans- 
verse eigenmode. This decomposition can be generalized 
for waveguides with confinement frequencies u>± that are 
slowly varying with the longitudinal coordinate x; in that 
case, the eigenstates Xn( r -iJ would vary with x as well, 
through their parametric dependence on uj± [371 ] - 

In general, the coherent source in Eq. ^ populates a 
finite linear combination of transverse eigenmodes in the 
waveguide, which would be coupled to each other through 
the interaction term. Clean bosonic scattering experi- 
ments, however, would ideally require a well-defined lon- 
gitudinal kinetic energy of the atoms, which means that 
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FIG. 1: (Color online) Possible waveguide geometry for a 
guided atom laser experiment with Bose-Einstein condensates 
in the presence of an optical scattering potential. The lower 
panel shows the effective potentials along the longitudinal 
coordinate x for the reservoir (V r ; green parabola) and the 
waveguide atoms (V g ; black horizontal curve). A radiofre- 
quency field with suitable frequency Lu T f (indicated by the ver- 
tical red line) can be used to resonantly couple atoms from 
the trapped condensate (whose chemical potential is denoted 
by the horizontal blue line) into the waveguide. 



the population ought to be restricted to one single trans- 
verse mode. Technically, this task could possibly be ac- 
complished by imposing a smooth barrier potential along 
the waveguide with a maximum height V" max that satis- 
fies fi — 2Hlu± < Vmax < /J, — fkuj_. In this case, atoms in 
the transverse ground mode (with the offset energy fkuj_ ) 
would have just enough kinetic energy to pass that bar- 
rier, while components associated with excited transverse 
modes (with offset energies > 2hui±) would be reflected. 

In the following, we assume that the population of the 
waveguide by the source can indeed be restricted to the 
transverse ground mode Xo( r ±) ex p(— muj±r' 2 L /2h). 
We then make the simplifying ansatz 



$ fl (r,t) = Xo(r±)ip(x,t) 



(8) 



which restricts the consideration to this ground mode 
and neglects the effect of (possibly virtual) transitions to 
excited transverse modes due to the interaction. ip{x, t) 
then satisfies the equation 

ih—${x,t) = H^ x) i)(x,t) + gft{x,t)${x,t)i;{x,t) 

+S(x,t)exp(-ifd/h) (9) 
with the one-dimensional single-particle Hamiltonian 

h 2 d 2 



H, 



(x) 



+ V(x) + huj±_ 



(10) 



2m dx 2 

and the effective one-dimensional interaction strength 
muj± 

g = ~2^ U a = 2nw ± a s (11 
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where a s is the s-wave scattering length for atoms in the 
state \g). The source amplitude in this one-dimensional 
equation is given by 

S(x,t) = J d 2 r ±X *(r ± )y r (r)K(t). (12) 

Assuming that the condensate wavefunction ty r (r) has no 
overlap with the region in which the scattering potential 
V(x) is defined, we make in the following the idealized 
ansatz of a strongly localized source, namely 

S(x,t) = S (t)8(x - x ) (13) 

where xq corresponds to the position around which the 
trap is centered. 

Our aim is to compute the time evolution of the con- 
densate wavefunction in the waveguide as well as the 
amount of depletion that is generated during this trans- 
port process. To this end, we follow the lines of the treat- 
ment in Refs. [12, [33| and make the Bogoliubov ansatz 
for the field operator tf}(x,t) in the transverse ground 
mode of the waveguide, which is decomposed into its or- 
der parameter (ip(x,t)} and a quantum fluctuation oper- 
ator Sip(x,t) according to 

$(x ) t) = $(x,t))+8H(x,t). (14) 

Using Eq. ([9]), it is then straightforward to de- 
rive an infinite set of coupled kinetic equations for 
the condensate wavefunction (tf>(x, i)) and the expec- 
tation values of products of the operators 8ij>(x,t) 
and their hermitean conjugates - - i.e., for the n- 
point functions (8ip(xi,t)8i/j(x2, t)), (8ij)'(xx, t)8ip(x2, t)), 
(54>(xi, t)Sip(x2, t)Sip(x3, t)), etc. The time evolution of 
(ip(x,t)), for instance, is described by the equation 

inl{^(x,t)) = H ( x) $(x, t)) + g\$(x, t))\ 2 (ip(x, t)) 

+2g(d^(x,t)5i>(x,t))(i>(x,t)) 
+g(8j>(x,t)5$(x,t))(ft{x,t)) 
+g{8${x,t)8i){x,i)8i){x,t)) 
+S(x,t)esq?(-ifit/K), (15) 

which contains couplings to n-point functions up to n — 
3. The equations for the higher-order cumulants, e.g. for 
(8ip(xi, t)8tp(x2, t)), are obtained by the time derivative 
of Sip(x, t) = tp(x, t) ~ (ip{x, t)) through the combination 
of Eqs. flSj) and (fl5|. 

In order to define a meaningful truncation scheme for 
this infinite hierarchy of equations, we make the stan- 
dard mean-field assumption that we have a rather small 



interaction strength g, corresponding to a small s-wave 
scattering length a s of the atoms in the waveguide, and a 
comparatively large source amplitude So, giving rise to a 
large longitudinal density no = (m/2u)\So/h\ 2 of atoms 
that are injected by the source. The mean-field regime 
could then be expressed by the inequalities no A 3> 1 and 
a s /\ <C 1 where A characterizes the relevant length scales 
of the waveguide system [such as the transverse oscilla- 
tor length a± = y/H/mu)±, the de Broglie wavelength of 
the incident matter-wave beam, or the length scales in- 
troduced by the scattering potential V(x)]. Formally, it 
would be defined by the simultaneous limits hq — ■» oo and 
a s — * with a s no remaining constant (38j |. 

This formal mean-field limit allows us now to classify 
the individual terms in Eq. (|15p and in the equations for 
the higher-order n-point functions according to different 
powers in the large parameter N = noA characterizing 
the number of atoms wi thin the length scale A. As we 
have | So | = hyj2una/m cx \/N, we can obviously infer 
from Eq. (fT5|) that the condensate wavefunction (ip(x, t)) 
scales as 0(N 1 ^ 2 ) for finite evolution times. The equa- 
tion for (8tl){xi,t)8ip(x2,t)) [see Eq. (|20)) below] is not 
directly connected to the coherent source, but contains 
the inhomogeneous term g8{x\ — X2){tp(xi,t)){if)(x2,t)) 
through which this two-point function becomes pop- 
ulated in the course of time evolution. Taking into 
account that g cx l/N in this formal mean-field 
limit and using (ip(x,t)) ~ 0(N 1 ^ 2 ), we then obtain 
(8ip(xi,t)8ip(x2,t)) ~ O(N ). Similarly one can show 
(8^(x 1 ,t)8^(x2,t)) ~ O(N ) as well, while three-point 
functions such as (8?p(xi,t)8'>p(x2,t)8'ip(x3,t)} would 
scale as ©(A^- 1 / 2 ). 

Following the notation introduced by Kohler and Bur- 
nett [33|], we define the condensate wavefunction 

iP{x,t) = {i>(x,t)}, (16) 

for atoms in the transverse ground mode of the waveg- 
uide, the pair function 

$(xi,x 2 ,t) = {8tp(x 2 ,t)8ip(xx,t)), (17) 

and the one-body density matrix of noncondensed atoms 

r(a;i,a:2,*) = (5ft(x 2 ,t)6$(xi,t)). (18) 

By taking into account all terms that scale at least as 
0(N^ 1 ), we obtain a closed system of coupled equations 
for ip(x, t), &(xi, X2,t), and T(x\, X2,t), namely 
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d 



ih—$(x 1 ,x 2 ,t) 



ifi-^T(xx,x 2 ,t) 
at 



H^ x) ip(x, t)+g \ip(x, t)\ 2 ip(x, t) + S(x, t) exp(-ifit/K) 
+2g^(x, t)T(x, x, t) + gip* (x, t)$(x, y, t) 



H, 



■2g\il>( Xl ,t)\' 



to) 



2g\^(x 2 ,t)\ 2 $(xx,x 2 ,t) 



+gS(xi - x 2 )ip(xx,t)ip(x 2 ,t) 
+gip 2 (x 2 ,t)T(xx,x 2 ,t) + gip 2 (x 1 ,t)T(x 2 ,x 1 ,t) 

+g5(x 1 - x 2 )$(x 1 ,x 2 ,t) + 2g [T(x 1 ,x 1 ,t) + Y{x 2 ,x 2 ,t)] $(xx,x 2 ,t) 
+gT(x 2 , x x , ,x 1 ,t)+gT(x 1 ,x 2 , t)$(x 2 ,x 2 ,t) 

H^ Xl) + 2g \ip(xx,t)\ 2 - Ht 2) - 2g \^(x 2 , t)\ 2 ] T{x x ,x 2 ,t) 

+gi> 2 {x 1 ,t)<P*(x 1 ,x 2 ,t) - g (*P*(x 2 ,t)) 2 ®(xx,x 2 , t) 

+2g \T(xx,Xi,t) - T(x 2 ,x 2 , £)] T(xx,x 2 ,t) 

+g$* (xx , x 2 , t)$(xx , xx , t) - 5$ (x x , x 2 , f )$* (x 2 ,x 2 ,t) 



0(N 1 ' 2 ) 
0(N- 1/2 



(19) 



[0(N )] 
[OiN- 1 )} (20) 

[0(N°)] 

[0(N~ 1 )] (21) 



r 



where consistently all couplings to higher-order n-point 
functions [which would scale as 0(iV~ 3 / 2 ) in Eq. (TIT))) and 
as 0(N- 2 ) in Eqs. (|2QI2ip ] are omitted. 

This set of kinetic equations is equivalent to the 
second-order Hartree-Fock-Bogoliubov approximation in 
modal form [35| . It is norm- conserving, i.e., for x 7^ xq 
the continuity equation 



l-n(x,t) + ^j(x,t) = 

is satisfied with 

n(x,t) = (x, f)V>(a:, f)) 

= \iP{x,t)\ 2 +T{x,x,t) 

the longitudinal particle density and 



(22) 



(23) 



h 

2mi 



r{x,t)^{x,t)-^{x,t)^-r(x,t) 



d 



d 



dxx dx 2 



T(xx,x 2 ) 



X1—X2—X 



(24) 



the atomic current along the waveguide. As for the den- 
sity n(x, i), the current j(x, t) consists of a coherent com- 
ponent associated with the condensate and an incoherent 
component that represents the noncondensed fraction, 
given by the second and third line of Eq. (|24|) . respec- 
tively. 

The transmission of a scattering state that is associ- 
ated with a given source amplitude So is obtained by 
dividing j(x, t) through the incident current ji. The lat- 
ter is defined as the current that would be obtained at 
the same source amplitude So if both the scattering po- 
tential V(x) and the interaction strength g were zero. We 



straightforwardly obtain for that case the asymptotic so- 
lution 



ip(x,t) = -^A/^exp 
h y 2/i 



^ (^/2m,Li\x - x \ - \xt 



which yields 



Ji 



(25) 
(26) 



We point out that the interaction strength g is assumed 
to vanish at a; = Xo, which avoids effects of nonlinear 
backaction with the source. In between the source and 
the scattering potential, g is adiabatically increased with 
1 up to a given maximal value, which remains then con- 
stant along the support of V(x). 

For the numerical treatment of the kinetic equations 
(|19H21[) , we employ a finite-difference representation of 
the one- and two-point functions t/>(x, t), $(xi, x 2 , t), and 
r(xi, x 2 , t) 26]. The numerical integration of the result- 
ing vectorial differential equations of the form iy = TLy is 
then performed by means of the Crank- Nicolson method 
[3^], which utilizes matrix- vector products as well as 
the solution of linear matrix-vector equations involving 
the effective Hamiltonian matrix TL. The split-operator 
technique is applied in the equations for $ and T in 
order to do separate propagation steps in the x\ and 
x 2 coordinates. To account for the nonlinearities, i.e. 
the dependence of H on y, each integration step is 
"corrected" by repeating it with H being evaluated at 
[y(t)+y {1) { t + St )}/ 2 where Y^U + St) denotes the first- 
order estimate for y(i + St) [26ll4fj|. 

In this grid representation, the Dirac delta functions 
appearing in the source term (|13p and in the interac- 
tion [cf. the second line of Eq. (|20p ] are naturally re- 
placed by Kronecker deltas with the appropriate prefac- 
tors. We verified that the regularization of the interaction 
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potential by a narrow Gaussian function does not signif- 
icantly change our numerical results. The amplitude of 
the source is adiabatically ramped from zero to So within 
< t < to according to 



S (t) 



7T t 

sin | — — 

2 t 



(27) 



To avoid nonlinear backaction of the condensate with the 
source, we also introduce a spatial ramping of the inter- 
action strength g according to 



1 



g(x) = — [1 + tanh (n(x — x%))] g 



(28) 



where n and x\ are chosen such that g(xo) < O.Olg at 
the position xq of the source and g(x) > 0.99g within 
the spatial region in which the scattering potential V(x) 
is finite. We verified that the presence of this ramping, 
which in practice could be induced by a spatially varying 
confinement frequency luj_ = to±(x), does not cause a 
significant additional amount of reflection. 

To avoid artificial backreflection of atoms from the 
boundaries of the numerical grid, we furthermore intro- 
duce complex absorbing potentials of the form V a f, s (x) = 
—ij(x) where 7(x) is vanishing throughout the scattering 
region and adiabatically ramped to a finite positive value 
in the vicinity of the boundaries [U \^ . As the single- 
particle Hamiltonian of the system thereby becomes non- 
hermitean, care is required for the inclusion of this imag- 
inary potential into the equation of the one-body density 
matrix T(xi,X2)', we obtain — «[ , y(a;) + j(y)]T(x, y, t) as 
an additional term in Eq. (f2"Tj) . 



III. TRANSPORT THROUGH DOUBLE 
BARRIER POTENTIALS 

Let us first investigate resonant transport of a Bose- 
Einstein condensate through a symmetric double barrier 
potential, which can be interpreted as a Fabry-Perot res- 
onator for coherent matter waves [l8l . [l9l |. Such an atomic 
quantum dot may exhibit long-lived quasi-bound eigen- 
states that allow for perfect transmission at energies in 
the immediate vicinity of the eigenenergies of such states, 
while off-resonant components of the incident wave will 
undergo reflection. In the context of guided atom lasers, 
such double barrier potentials may therefore serve as 
tools for energetic purification of the incident beam, in 
the case that the source is, in practice, not yielding a 
perfectly monochromatic matter wave. 



In accordance with our previous studies [24j |. we con- 
sider a double barrier potential of the form 



V(x) = V Q 



-(x- Xl -L/2) 2 /a 2 , e -(x- Xl +L/2) 2 /a 2 



(29) 



with x\ the central position between the two barriers, 
L the distance between the maxima, and a the width 
of the barriers. The Bose-Einstein condensate is as- 
sumed to consist of 87 Rb atoms with the scattering length 




1.4 




1.15 



FIG. 2: (Color online) Time-dependent sweep across a res- 
onance of the double barrier potential ([29} [which is dis- 
played in the lower panel of Fig. (J3J]. The upper panel 
shows the nonlinear transmission spectrum at the incident 
current ji = 1.6u>± obtained from the integration of the 
Gross-Pitaevskii equation (|30jl with the interaction strength 
g ~ 0.034:hu!±a±. Dashed lines mark the resonance branches 
that are not directly accessible by means of a straightforward 
injection of matter waves at fixed chemical potential. The 
dotted line represents the linear transmission spectrum at 
g = 0. The middle and lower panels show the time-dependent 
transmission that is obtained by sweeping the chemical po- 
tential according to Eq. (|3ip during the propagation process. 
This sweep allows one to populate the upper branch of the 
resonance peak that corresponds to the 5th excited quasi- 
bound state within the cavity. Apart from a slight shift of 
the distorted resonance peak (due to a renormalization of the 
effective interaction strength), the calculation based on the 
dynamical Hartree-Fock-Bogoliubov (HFB) equations (|19H21[) 
(lower panel) reproduces quite well the Gross-Pitaevskii (GP) 
calculation (middle panel). 



a sc = 5.77nm and to propagate in a waveguide with the 
transverse confinement frequency ujj_ = 2ir x 10 s . In 
terms of uj±_ and the transverse oscillator length a±_ = 
^/h/moj± ~ 0.34/xm, we therefore have g ~ 0.03Ahoj±a±. 
As in Ref. [2H, we choose Vb = huj± and L = 10a = 5^m 
~ 14.7aj_ for the parameters of the potential |2"9"1) . 

The main characteristics of this transport process can 
be understood on the level of the inhomogeneous Gross- 
Pitaevskii equation 



. d . . 
9i 



= #< X V(M) + <?I#M)| 2 #M) 



-S(x, t) exp(— i[it/h) 



(30) 



which corresponds to the lowest-order truncation of the 
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system of equations (fT9H2H . The upper panel of Fig. [2] 
displays the transmission spectrum through this double 
barrier potential as a function of the condensate's chem- 
ical potential fi, in the absence and presence of the in- 
teraction between the atoms (dotted and solid/dashed 
lines, respectively), which is computed by means of sta- 
tionary scattering solutions of the Gross-Pitaevskii equa- 
tion that exhibit the incident current ji = 1.6ojj_' For 
a noninteracting condensate, a sequence of Breit-Wigner 
peaks would be obtained in the transmission, whose po- 
sitions and widths correspond to the energies and decay 
rates, respectively, of the quasi-bound states within the 
resonator. These transmission peaks exist also in the 
presence of interaction, but become strongly distorted 
towards higher chemical potentials (for repulsive inter- 
action) due to the nonlinearity in the Gross-Pitaevskii 
equation [U, [26| . A multivalued transmission spectrum 
with several bistable branches is thereby obtained, which 
is typical for nonlinear transmission problems and arises 
also in nonlinear optics [3l[ as well as in the electronic 
transport through quantum wells [H, 0, [4f| . 

It is easy to see that a direct, straightforward injec- 
tion of matter waves onto this atomic Fabry-Perot res- 
onator at a given chemical potential fi would always lead 
to the population of the lowest branch in the transmission 
spectrum [24l [26| . Hence, resonant transmission through 
long-lived quasi-bound states with narrow widths can- 
not be achieved in this way, and the atomic quantum 
dot qualitatively acts like a simple potential barrier, giv- 
ing rise to nearly perfect (classical-like) reflection of the 
condensate for chemical potentials well below the barrier 
height Vq. Indeed, if the chemical potential of the in- 
cident beam is chosen to match the level of an internal 
(noninteracting) quasi-bound state, a finite population 
of this state will shift the associated level towards higher 
energies, and the matter-wave beam will no longer be on 
resonance. On the other hand, injecting the condensate 
at the appropriate chemical potential that corresponds 
to the position of the distorted resonance peak would not 
yield any population in the internal quasi-bound state in 
the first place if the resonator was initially empty. 

In Ref. [lij we proposed a control scheme that allows 
one to overcome this limitation and to achieve resonant 
transport of the condensate in the presence of finite in- 
teraction. The basic idea is that the chemical potential 
of the condensate (or, equivalently, the offset potential 
of the waveguide) needs to be adiabatically varied dur- 
ing the propagation process, in order to follow the upper 
branch of one of the distorted resonance peaks. The fea- 
sibility of this scheme was demonstrated through numeri- 
cal simulations on the basis of the time-dependent Gross- 
Pitaevskii equation. Starting at a chemical potential on 
the left shoulder of the resonance peak corresponding to 
the 5th excited quasibound state (which is magnified in 
the middle and lower panels of Fig. [2]), and increasing 
the chemical potential up to the value that equals the 
peak's top position, we obtained nearly perfect transmis- 
sion on finite time scales, limited only by the dynamical 



instability of the nonlinear scattering state at resonance 

The middle and lower panels of Fig. [5] display the 
outcome of a similar time-dependent transport process, 
which was calculated here with the Gross-Pitaevskii 
equation (middle panel) as well as through the the nu- 
merical integration of the dynamical Hartree-Fock Bo- 
goliubov equations (fT9H2"T)) . In these calculations, we 
first populated a scattering state at the chemical poten- 
tial fio — 0.97fnu±, by ramping up the source amplitude 
within the propagation time to — lO 3 ^ 1 , and then sub- 
jected it to a linear sweep of the chemical potential ac- 
cording to 

H(t) =/*> + («- to)£ (31) 

with (i = 5x 10 -5 fia>^. Both calculations lead to a similar 
scenario, namely the adiabatic population of the upper 
branch of the resonance peak corresponding to the 5th 
excited quasi-bound state in the resonator, and the sub- 
sequent decay to a low-transmission state at propagation 
times t at which the corresponding chemical potential 
/i(t) exceeds the position of the resonance peak. A slight 
displacement in the peak position and in the slope of the 
resonance branch is observed in the HFB calculations 
when being compared with the nonlinear transmission 
curve obtained from the Gross-Pitaevskii equation. We 
attribute this to a renormalization of the effective non- 
linearity strength in Eqs. (fl9H2T|) due to the coupling of 
the condensate wavefunction to the two-point functions 
$(xi,X2,t) and r(xi,X2,i). 

Fig. [3] displays the condensate density (n c ) as well as 
the density of noncondensed atoms (n nc ) of the instan- 
taneous scattering state at the propagation time t with 
fi(t) — 1.1 at which an overall transmission of T ~ 0.85 is 
obtained. We clearly see a finite cloud of noncondensed 
atoms within the resonator, which is obviously induced 
by the large local density fluctuations of the condensate. 
The noncondensed denisty is found to slowly increase 
with time during the sweep of the chemical potential, but 
remains rather low compared to the condensate density 
on the time scales that are considered here. 

Our study indicates that the control scheme pro- 
posed in Ref. [HI is working not only within the Gross- 
Pitaevskii description of the transport process, but re- 
mains feasible also from a more microscopic point of view 
taking into account the noncondensed cloud that is gen- 
erated within the scattering region. We remark, how- 
ever, that this conclusion strongly depends on the specific 
parameters at hand. Indeed, if we increase the micro- 
scopic interaction strength of the atoms by a factor two 
(which could possibly be achieved through a Feshbach 
resonance) while keeping the other parameters of the sys- 
tem fixed, we find rather strong depletion of the con- 
densate within the resonator during the time-dependent 
sweep of the chemical potential across the resonance. In 
this situation, the truncation scheme leading to the equa- 
tions (TlT)H2"Tj) would not be justified any longer, and we 
expect that at that particular high value of the interac- 
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FIG. 3: (Color online) Density of the condensate wavefunc- 
tion (upper panel) and of the noncondensed atoms (middle 
panel, also shown as red dashed line in the upper panel) near 
a resonance of the double barrier potential (|29[) . The plot 
represents a snapshot of the time-dependent sweep across the 
resonance shown in Fig. [2] at the evolution time at which the 
chemical potential equals /i = 1.1, corresponding to a trans- 
mission of T ~ 0.85. The lower panel shows the double barrier 
potential. 



tion strength the overall coherence of the atomic cloud 
cannot be maintained and the condensate would become 
destroyed during the transport process. 



IV. TRANSPORT THROUGH DISORDER 
POTENTIALS 

Let us now consider the transport of Bose-Einstein con- 
densates through disorder potentials of finite range, con- 
sisting of a more or less random sequence of several bar- 
riers and wells with different heights and depths. Such 
disorder potentials received special attention in the past 
years due to the possibility of reali zing Anderson local- 
ization with ultracold atoms @, @, Hf |H EE El H3] • 
This task was ultimately achieved in two independent 
experiments [H, Q studying the expansion of condensates 
within optical speckle fields [f| as well as within bichro- 
matic optical lattices |], in a parameter regime where 
effects of the atom-atom interaction on the expansion 
process could be neglected. Indeed, various theoreti- 
cal studies revealed that the presence of a finite non- 
linearity in the wave equation that describes the expan- 
sion process, which would correspond to a finite inter- 
action strength within the mean-field description of the 
condensate, would lead to a significant deviation from 
the scenario of Anderson localization, involving indefi- 
nite spreading of the tails of the condensate wavefunction 
instead of exponential localization [5l|, [H, HH ■ 



The effect of a finite interaction is predicted to be even 
more dramatic in scattering processes through disorder 
potentials. Here, Anderson localization would manifest 
in an exponential descrease of the average transmission 
T with the spatial length L of the disorder region (or, 
more precisely, in a linear increase of the average (— InT) 
with L), which is induced by the destructive interference 
between multiply reflected components of the scattering 
wavefunction [551 ]. In the presence of a finite nonlinear- 
ity, however, a crossover to an algebraic decrease of the 
average transmission according to (T) ~ (L + i ) _1 is 
encountered (25| . which would essentially correspond to 
Ohm's law for incoherent transport. We find that this 
crossover is directly connected with the appearance of 
permanently time- dependent scattering of the condensate 
[25j, |56j. That is, the injection process of the condensate 
according to the working principle of the guided atom 
laser does not lead to the population of a quasi-stationary 
scattering state, but induces time-dependent fluctuations 
of the condensate density beyond a certain critical value 
of the source amplitude. 

The main finding of our investigation of this prob- 
lem with the kinetic Hartree-Fock-Bogoliubov equations 
([T5HH]) is that the appearance of permanently time- 
dependent scattering within the Gross-Pitaevskii descrip- 
tion of the transport process involves strong depletion 
and the production of a comparatively large cloud of 
noncondensed atoms. This indicates that the overall co- 
herence characterizing the condensate becomes destroyed 
on a microscopic level, which would indeed be consis- 
tent with an Ohm-like scenario of incoherent transport 
through the disorder region. We exemplify this behaviour 
with one particular, relatively short disorder potential, 
displayed in the lower panels of Figs. [4] and [5l For 
the chemical potential of the incident beam, we consider 
fi = 1.0hu± which is of the same order as the barrier 
heights of the disorder potential. As in Section IIII1 we 
choose u>± = 2ir x 10 3 s _1 for the transverse confinement 
frequency of the waveguide and g = 0.034huj±a± for the 
effective one-dimensional interaction strength. 

Let us first consider the transport process across this 
disorder potential at the incident current ji = 1.5wj_ 
for which a Gross-Pitaevskii calculation yields quasi- 
stationary scattering. At this relatively low value of the 
incident current, we find nearly stationary dynamics on 
the level of the kinetic Hartree-Fock-Bogoliubov equa- 
tions (|19r[2Tj) as well, where the density of the scattering 
state is, as shown in Fig. [5J dominated by the coher- 
ent component associated with the condensate. As for 
the double barrier potential discussed in Section IIII1 a 
finite, slowly increasing amount of depletion is generated 
within the scattering region, with noncondensed atoms 
preferrably accumulating around the combined minima 
of the disorder potential and the condensate density. On 
time scales of the order of lOOOo^J 1 corresponding to 160 
ms, however, the total density of noncondensed atoms is 
still rather low compared to the condensate, which indi- 
cates that the coherence of the atomic beam should be 





FIG. 4: (Color online) Density of the condensate wavefunc- 
tion (upper panel) and of the noncondensed atoms (middle 
panel, also shown as red dashed line in the upper panel) 
during the time-dependent propagation process of the con- 
densate through a disorder potential, at the evolution time 
t = 1000/cjx (corresponding to 160 ms). The incident current 
ji — 1.5o>± was chosen such that a Gross-Pitaevskii calcula- 
tion of the transport process yields a quasi-stationary scat- 
tering state. During the propagation process, noncondensed 
atoms are found to slowly accumulate preferrably around the 
minima of the disorder potential, which is shown in the lower 
panel. The total depletion, however, is still relatively low at 
t = 1000/wx- 



well preserved also on the microscopic level. 

The situation is rather different for a high incident 
current such as ji — 4.775wx, for which the simula- 
tion of the scattering process on the basis of the Gross- 
Pitaevskii equation yields permanently time-dependent 
fluctuations. Such time-dependent behaviour is also 
found on the level of the kinetic Hartree-Fock-Bogoliubov 
equations (fl9H2T|) . together with a rather large amount 
of depletion, as shown by the snapshot displayed in 
Fig. [5l In that case, the truncation scheme lead- 
ing to the closed system of equations (fl9rl21"j) can no 
longer be justified, and higher-order cumulants, such as 
(8i/j(xi,t)5ij)(x2,t)5il)(x3,t)) for instance, would have to 
be taken into account as well for an accurate descrip- 
tion of the transport process. We conjecture that in this 
situation the condensate fraction becomes destroyed on 
a microscopic level, and the transport process becomes 
incoherent. 

Fig. [6] provides a comparison of the time-dependent 
currents of the condensate and the noncondensed atoms 
[cf. Eq. (pMj) ] for the cases of low and high intensity of the 
incident matter-wave beam. While a stable stationary 
flow, dominated by the coherent component of the con- 
densate, is encountered in the presence of a low incident 
current, time-dependent dynamics arises already during 



FIG. 5: (Color online) Same as Fig. [3] for a higher incident 
current ji — 4.775lu± for which the Gross-Pitaevskii calcula- 
tion yields permanently time-dependent scattering. The fig- 
ure shows a snapshot at the evolution time t — 500 / u>± . Here 
the density of noncondensed atoms n nc is of the same order as 
the condensate density, which indicates strong depletion and 
the destruction of the condensate on the microscopic level. 




500 



FIG. 6: (Color online) Total current of the condensate and 
the noncondensed atoms (blue and red lines, respectively) as 
a function of the propagation time, for the incident currents 
ji = 1.5lo± (upper panel) and ji = 4. 775a; x (lower panel), in 
the presence of the disorder potential that is shown in Fig. [4] 
The plots include the time interval of the initial ramping pro- 
cess of the source amplitude, which lasts until t ~ 300/a;x. 
In the case of a rather low incident current (upper panel) 
for which a Gross-Pitaevskii calculation would yield quasi- 
stationary transport, depletion remains rather low on finite 
time scales. At higher incident currents for which perma- 
nently time-dependent scattering manifests on the Gross- 
Pitaevskii level, we find strongly enhanced depletion and a 
current of noncondensed atoms that is comparable to the cur- 
rent of the condensate fraction. This indicates, on a micro- 
scopic level, that the condensate becomes destroyed during 
this time-dependent scattering process. 
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FIG. 7: Total atomic current across the disorder potential as 
a function of time, calculated for the critical incident current 
ji = 2.387ojx at which time-dependent scattering begins to 
appear on the level of the Gross-Pitaevkii equation (upper 
panel). A simulation of this scattering process with the ki- 
netic Hartree-Fock-Bogoliubov equations (lower panel), how- 
ever, yields quasi-stationary behaviour on a time scale of the 
order of 1000a/[ , with a reasonably low amount of depletion 
(between 15 and 20 percent of the total density), while dy- 
namical instability and strong depletion sets in at ui±t ~ 1500. 



the ramping process of the source for ji — 4.775cl>j_ and 
induces in turn rather strong depletion from u±t ~ 350 
on. Obviously, our calculation is, in this latter case, not 
trustable any longer beyond that evolution time as far as 
the accurate prediction of the time-dependent densities 
and currents of the condensate and the noncondensed 
fractions is concerned. 

The scenario discussed here is systematically found 
also for a number of other disorder potentials of similar 
length as the one shown in Figs. 2] and Stationary scat- 
tering on the Gross-Pitaevskii level implies, in general, a 
nearly stationary and fairly coherent flux on the level of 
the kinetic Hartree-Fock-Bogoliubov equations, whereas 
time-dependent scattering within the Gross-Pitaevskii 
equation leads to strong depiction. We point out, how- 
ever, that there is no rigorous one-to-one relation be- 
tween these two phenomena. Indeed, exceptions from 
this general rule particularly arise if the incident current 

(c) 

of the condensate is close to the critical value j\ beyond 
which permanently time-dependent scattering occurs on 
the level of the Gross-Pitaevskii equation. In this case, 
a small modification of the effective nonlinearity, caused 
by the coupling to the noncondensed cloud, can have a 
drastic effect and may destabilize — or, in some cases, 
also re-stabilize — the flow of the condensate. An exam- 
ple is shown in Fig. where the time evolution of the 
current across the disorder potential under consideration 
is shown for the critical incident current ji = 2.387oj^ 
beyond which the flux becomes unstable on the Gross- 
Pitaevskii level. In contrast to the Gross-Pitaevskii cal- 
culation, the simulation of this transport process with the 
kinetic Hartree-Fock-Bogoliubov equations does not give 
rise to time-dependent fluctuations of the current on time 
scales of the order of lOOO^ 1 (which would correspond 
to 160 ms) and leads to a stable flux with a reasonably 



low amount of depletion, until destabilization sets in later 
on at u>±t ~ 1500. 



V. CONCLUSION 

In this work, we investigated the transport of guided 
bosonic matter waves through one-dimensional scattering 
potentials by means of the inhomogeneous Hartree-Fock- 
Bogoliubov equations, which simulate the coherent injec- 
tion of those matter waves from a trapped Bose-Einstein 
condensate [l2j]. A straightforward numerical integra- 
tion of these Hartree-Fock-Bogoliubov equations yields 
nearly stationary, mean-field-like scattering dynamics of 
the bosonic flow on reasonably long time scales, provided 
the effective one-dimensional interaction strength is not 
too large and the integration of the corresponding Gross- 
Pitaevskii equation leads to quasistationary scattering. 
In that case, a comparatively low amount of depletion is 
obtained, with noncondensed atoms preferrably accumu- 
lating around the local minima of the scattering poten- 
tial. If, on the other hand, the nonlinear scattering dy- 
namics of the condensate is explicitly time-dependent on 
the level of the Gross-Pitaevskii equation, we encounter 
strong depletion and a rather large density of noncon- 
densed atoms within the Hartree-Fock-Bogoliubov de- 
scription. This indicates that the condensate becomes 
destroyed on a microscopic many-body description of the 
transport process, and that the atoms are individually 
propagating across the scattering potential. 

It would be interesting, though beyond the scope of 
this work, to compare our findings with the truncated 
Wigner approach |57l . [H, [5i|, which represents an al- 
ternative method to describe dynamical transport pro- 
cesses of Bose-Einstein condensates beyond the Gross- 
Pitaevskii description. In this approach, the microscopic 
quantum dynamics of the atomic cloud is described by a 
set of classical fields which sample the distribution of the 
initial many-body state describing the Bose gas under 
consideration, and which evolve according to the time- 
dependent Gross-Pitaevskii equation. In our situation, 
this initial state would correspond to an almost per- 
fect Bose-Einstein condensate of the "reservoir" atoms 
that are confined in the three-dimensional magnetic trap. 
The classical fields that would have to be used for the 
truncated Wigner approach would therefore be rather 
similar at the initial stage of the matter-wave injection 
process, and deviations between them, which ultimately 
give rise to a loss of coherence, should mainly arise in 
the presence of strongly time-dependent scattering where 
the occurrence of classical chaos in the Gross-Pitaevskii 
equation might induce an exponential sensitivity of the 
time evolution with respect to initial conditions. This 
would be consistent with our observation that depletion 
is particularly strong in the regime of permamently time- 
dependent scattering. 

An interesting study in this context was recently re- 
ported in Ref . [60j , where the quantitative amount of de- 
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pletion was computed for a homogeneously moving three- 
dimensional condensate (initially defined as a plane wave 
in a box with periodic boundaries) in the presence of 
an obstacle or a disorder potential. Using the truncated 
Wigner method, the authors showed that nonstationary, 
turbulent flow of the condensate leads to rather rapid 
depletion, on time scales that seem to be roughly compa- 
rable with our investigations (olTj . Good agreement be- 
tween the time-dependent Hartree-Fock-Bogoliubov ap- 
proach and the truncated Wigner approach was, more- 
over, reported in other studies [62| that were focusing on 
collapsing Bose-Einstein condensates. 

We finally point out again that several idealizations 
and approximations were employed for the desciption of 
the bosonic transport process. Specifically, we assumed 
a perfect condensate in the reservoir, which is not sig- 
nificantly depopulated during the outcoupling process, 
we neglected the (virtual) population of transversally ex- 
cited modes in the waveguide, and we excluded nonlin- 
ear back-action with the source assuming that the inter- 
action strength of the waveguide atoms is negligible in 
the vicinity of the magnetic trap. These simplifications 



would have to be reconsidered for the simulation of a re- 
alistic matter-wave propagation experiment based on an 
atom laser. We do, however, not expect that the general 
findings obtained in this study would qualitatively alter 
if these complications are properly taken into account. 
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